Phase diagram and density large deviations of a nonconserving ABC model 
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The effect of particle-nonconserving processes on the steady state of driven diffusive systems is 
studied within the context of a generalized ABC model. It is shown that in the limit of slow 
nonconserving processes, the large deviation function of the overall particle density can be computed 
by making use of the steady state density profile of the conserving model. In this limit one can define 
a chemical potential and identify first order transitions via Maxwell's construction, similarly to what 
is done in equilibrium systems. This method may be applied to other driven models subjected to 
slow nonconserving dynamics. 
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Driven diffusive systems have been at the focus of ex- 
tensive theoretical and experimental studies in recent 
years. Many studies have been devoted to systems of 
particles evolving under local biased exchange processes 
in the bulk which may or may not be coupled to non- 
conserving reservoirs at the boundaries. Examples in- 
clude the asymmetric exclusion process (ASEP), the 
zero-range process (ZRP) and many others [1- 4]. Bulk- 
nonconserving processes are rather common in many 
physical systems such as molecular motors, traffic flow 
problems with road intersections, chemical reactions in 
an open environment and others (see for example Q)- 
The corresponding driven models of such systems are of- 
ten less tractable analytically [5-12]. Their steady-state 
properties are usually very different from those of the 
corresponding conserving models. 

In this Letter we study the effect of bulk nonconserv- 
ing processes by considering the limit where they occur 
on a much longer time scale than that of the conserving 
dynamics. This separation of time scales allows the sys- 
tem to relax to its conserving steady state with a fixed 
number of particles in between nonconserving dynami- 
cal moves. As a result, the steady state of the system 
can be expressed in terms of an appropriate 'ensemble' 
of the conserving steady states which may be regarded as 
a kind of nonequilibrium grand-canonical ensemble. One 
may thus use the steady state properties of the conserv- 
ing system to analyze the corresponding properties of the 
nonconserving one. 

This approach is demonstrated on the ABC model 
il,[ii]. This is a one dimensional three species driven 
exclusion model which exhibits a phase separation tran- 
sition, and which has previously been generalized to in- 
clude particle-nonconserving processes. By considering 
the limit of slow nonconserving dynamics and following 
the approach described above, we compute explicitly the 
large deviation function (LDF) of the overall particle den- 
sity, despite the fact that the LDF of the density profile 
is not known. The LDF yields a definition for a 'chemi- 
cal potential' which unlike the equilibrium one, depends 
on the details of the nonconserving dynamical process. 



Based on this derivation, we draw the exact phase dia- 
gram of the model. The first order transition line exhib- 
ited by the model may be computed via Maxwell's con- 
struction, even though detailed balance is not obeyed. As 
discussed at the end of the Letter, the method presented 
below is readily applicable to other driven models that 
are coupled slowly to an external reservoir. 

The ABC model is defined on a one-dimensional peri- 
odic lattice of length L, where each site is occupied by 
one of the three species of particles, labeled A, B and 
C. The model evolves by random sequential updates 
whereby particles on neighbouring sites are exchanged 
with the following rates. 



AB ^ BA 
1 



BC ^ CB 

1 



CA ^ AC. (1) 
1 



For g = 1, the model relaxes to an equilibrium steady 
state where the particles are homogeneously distributed. 
For any finite value of q 7^ 1, the model exhibits phase 
separation into three domains in the limit oi L 00. 
Generically, the model does not obey detailed balance 
and it relaxes to a nonequilibrium steady-state. A spe- 
cial feature of the ABC model is that in the case where 
the number of particles of the three species are equal, 
Na = Nb = Nc, the dynamics obeys detailed balance 
with respect to an effective Hamiltonian with long-range 
interactions. 

The model is often studied in the limit of weak asym- 
metry where q approaches 1 in the thermodynamic limit 
as q — exp(— [15]. This model exhibits a phase tran- 
sition at some value of f3 between a homogenous phase 
and an ordered phase with three macroscopic domains, 
each predominantly occupied by one of the species. The 
transition point and its order depend on the values of 
Na,Nb and Nc In the equal-densities case the 

parameter /3 plays the role of the inverse temperature. 
This phase transition has also been studied by consider- 
ing various generalization of the ABC model such as in- 
terval boundary conditions [13] , species-dependent q [13] 
and particle-nonconserving dynamics 

In the following we analyze the generalized weakly- 
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asymmetric ABC model on a ring with particle- 
nonconserving dynamics [l^ [2l[ . In this model sites can 
be occupied by inert vacancies, denoted by 0, whose dy- 
namics is defined as 

AO^OA, BO^OB, CO^OC. (2) 
111 

The total number of particles, N ^ Na + Nb + Nc < L, 
fluctuates through evaporation and deposition of triplets 
of neighbouring particles given by 



ABC 



p e 



000. 



(3) 



Here p is the overall rate of the nonconserving dynam- 
ics and fi is a parameter which as shown below, can be 
regarded as a chemical potential. This type of noncon- 
serving process has been chosen because it maintains de- 
tailed balance with respect to an effective Hamiltonian if 
initially the densities are equal, Na = Nb = Nc- 

In the equal-densities case, the equilibrium canonical 
and grand-canonical ensembles of this model correspond 
respectively to a conserving model, defined by rules ([!]) 
and ([2]) and a nonconserving model, which includes also 
rule ([3]). It has been found that for this case both mod- 
els exhibit the same second order transition line, which 
turns into a first order line at a tricritical point only in 
the nonconserving model [2^ 21 1. Such inequivalence be- 



tween the phase diagrams of the two ensembles is often 
observed in long-range interacting systems (for recent re- 
views see [23, [23 ) ■ 

In this Letter we study the phase diagram of the gen- 
eralized weakly-asymmetric ABC model with nonequal 
densities. In this case no free energy exists. This study 
can be accomplished by analyzing the steady state of its 
hydrodynamic equations, given by 



dtPa = I3L [pa {Pa+l - Pa+2) 

+ p{pl-e-^^^PAPBPc) 



L-'dip 



(4) 



Here pa {x) is the coarse-grained density profile whose 
index a denotes the species and runs cyclicly over A , B 
and C. We denote the average density of each species 
by ra = Na/ L — dxpa (x) and the overall density by 
r ^ rA + rB +rc- In the original ABC model (r = l,p = 
0); Eq. (|4]) has been shown to be exact in the limit of 
L 00 for equal densities (lil . [l7| . and has been argued 
to remain so even for arbitrary average densities [17,]. 

The conserving steady state (r < l,p = 0) of Eq. (jU, 
denoted by ix,r), can be readily extracted from the 
known steady state of the original ABC model, p* (x, 1), 
via a scaling transformation discussed below. In the non- 
conserving model, we are able to derive the steady-state 
of Eq. Q in the limit of slow nonconserving dynamics. 



where the p-dependent term becomes subdominant. Dy- 
namically this limit implies that on time scales of or- 
der the nonconserving model relaxes to the conserv- 
ing steady state with a fixed overall particle density, r, 
whereas on longer time scales of order L^, r fluctuates 
around its steady-state value. 

This separation of time scales suggests that Eq. (|4]) 
remains valid within the limit considered in Eq. ([5]) and 
that the steady state measure of the nonconserving model 
can be written in the limit of large L as 



(6) 



Here = {Q} denotes a microstate of the model with 
Ci = A,B,C for i G [1, i] and Pc (C ; N) is the conserving 
steady-state measure. Although Pc (C ; N) is not known, 
the knowledge of its extremizing profile in the hydrody- 
namic limit, /9* (x,r), is sufficient for deriving the prob- 
ability density of N in the nonconserving model, P{N). 

We now derive P{N) by writing its master equation, 
which evolves on the slow time scale of the bath by 

dtPiN) = P{N')Q{N' -> N)-P{N)Q{N -> iV')- 

N'=N±3 

(7) 

The average transition rates, Q, are given by 



Q{N^N') 



Pc{C-:N) W{C,N\C\N'), (8) 



E 



where N' — N zL 3 and W is the transition rate between 
two microstates. The latter is computed by counting 
the number of ABC and 000 triplets in each microstate, 
denoted by uabc (C) Eind nooo (C): respectively. In the 
limit of L — > (X) the evaporation rate can be computed 
using a saddle-point approximation, yielding 

Q{N^N-3) =pe-^^^ E ' ^) "^^BC (C) 

C 

~ p e-^f'i'L ( dxp\ (x, r) (x, r) pi (x, r) . (9) 



7>2, 



(5) 



The deposition rate is obtained similarly as 

Q{N^ N + 3)^pL f dx[p* (x, r)f , (10) 

JQ 

which can be simplified by noting that the inert vacancies 
have a fiat steady-state profile, Pq (x, r) = 1 — r. 

Equation ([7]) corresponds to a one-dimensional ran- 
dom walk in N in the presence of a local potential whose 
steady state is 

where N^i^ = N — Smiuc (Na) is the minimal number 
of particles. The LDF of the overall density, r, is 

Gip, r) ^-pr+ I dr'p (/) , 



(12) 
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Figure 1. The large deviation function of the overaU density, 
r, for P = 50, ta = tb = r/3 - 0.025 and fi = -0.053. For 
these parameters the LDF has two local minima, at r^am a-^d 
r*^.ci, corresponding to the homogenous and ordered phases, 
respectively. The nonconserving model undergoes a first order 
phase transition when G{^,r^^^) = G{fj,,r*^^). 



where tq is an arbitrary parameter and 



MO = [log 



dxp^p^'sp}, -3\ogil-r) . (13) 



The LDF, which is in fact proportional to the potential 
felt by the random walker, is plotted in Fig. [T] for some 
point in parameter-space. In the large L limit, the aver- 
age value of r is given by the global minimum of G (/i, r) 
which corresponds to dG/dr = jj, — p[r) ~ Q. 

We now discuss briefly the role of /i(r) as the chemi- 
cal potential of the conserving model. This will enable 
us to compare the phase diagrams of the two models in 
the (/3, /i)-plane. In the absence of a Hamiltonian which 
states the energy cost of adding or removing particles, the 
chemical potential of a conserving nonequilibrium system 
can be measured by coupling it to a microscopic gauge, 
customarily defined using Creutz method (2j|. Here, this 
is done by allowing triplets of ABC particles to evaporate 
from the lattice into a demon with the slow rate p defined 
in Eq. If the demon contains triplets, it may depose 
them back into the lattice at the same rate p. Following 
the lines of the derivation above, the probability density 
of the number particles in the demon, Nd > 0, can be 
shown to obey in the limit L — > oo. 



P(7V<i + 3) Jodxp*^p*sPh 



^ g3/3M('-) 



(14) 



and hence P {Nd) oc exp {/Sp. (r) Nd) UM- The function 
p (r) is therefore the chemical potential of the conserving 
model, as measured by the demon. In contrast to equilib- 
rium, here p{r) is derived from W{C, N\Q' , N') and thus 
depends on the choice of nonconserving dynamics. 

We now proceed to compute the phase diagrams of the 
model under conserving dynamics and slow nonconserv- 
ing dynamics. The conserving phase diagram (p = 0, r < 
1) can be derived from that of the original ABC model 
(p = 0, r = 1) using a mapping where the vacant sites are 
removed from each microstate of the conserving model. 
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Figure 2. The [1/ P, p) phase diagrams of the conserving (a) 
and nonconserving (b) models for two equal densities with 
A = 0.025. In the conserving model p is computed from Eq. 
(|13p . The thick and thin solid lines represent the first and 
second order phase transitions, respectively. They join at the 
conserving (•) and nonconserving (★) tricritical points. The 
dashed lines denote the stability limits of the two phases. The 
inset in (a) depicts schematically the area near the conserving 
tricritical point. The conserving tricritical point is irrelevant 
in (b), as it is located within the ordered phase of the non- 
conserving model. 



The master equation of the resulting system corresponds 
to that of the original ABC model with N sites. By ob- 
serving that q — exp{f3/L) — exp{/3r/N) = exp{(3' /N) 
we conclude that the iV-size system has an effective bias 
of (3' = /3r. Similarly the average densities of the TV- 
size system can be shown to be given by = Va/r for 
a = A, B,C and rg = 0. The steady-state profile of the 
conserving model can thus be expressed as 

P*a {x,l3,ra,r) = rp*^ {x,/3r,ra/r, 1) , (15) 

where p* in the r.h.s. has been derived for arbitrary Va 
in [16]. Equation (fT5|)_maps the phase diagram of the 
original ABC model [ISL llq onto the conserving model 
(r < l,p = 0). The resulting conserving phase diagram 
consists of a second order transition line at 



13 = 27r%/3/\/r2 - 36 A2, 



(16) 



r/3) is a measure for the 
The transition becomes 
r > 2(r\+r%+ r^) . The 



where A2 = i Ea=A,s,c " 
deviation from equal densities 
first order for (r^ + '^b + 
phase diagram is shown in Fig. [5^ for two equal densities, 

rA=rB^ r/3 - A, rc = r/3 + 2A. (17) 

Since the LDF of pa (x) in not known, one cannot com- 
pute the first order transition line of the conserving 
model. It is possible, however, to draw the stability limits 
(dashed lines) in between which both phases are stable 
with respect to small perturbations. 

The phase diagram of the nonconserving model can 
be derived by studying the extrema of G{p,r), given 
by the equation p = p{r). The function p{r), defined 
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Figure 3. The fi{r) curve for two equal densities with 
A = 0.025 and two values /3. The hydrodynamic solutions 
of the homogenous and ordered phases (solid lines) meet at 
the critical point (+). They are shown in comparison with 
the results of the conserving simulation (•) and nonconserv- 
ing simulation with homogenous (Kl) and fully ordered (A) 
initial states, performed with L = 2400 and p = 0.0001. The 
dashed line in (b) denotes the nonconserving first order tran- 
sition point. 



in Eq. is plotted for the two equal densities case 

in Fig. [21 At low values of (3 (Fig. [5^) there is a one- 
to-one correspondence between /i and r. The conserving 
and nonconserving models therefore behave similarly and 
display a second order phase transition. At high values 
of /3 (Fig. [3)3) we observe a region of fi where G (p, r) 
has three extrema. The intermediate density extremum 
has negative compressibility and corresponds to a maxi- 
mum of G (/i, r). It is thus stable only in the conserving 
model, while the nonconserving model undergoes a first 
order transition according to the Maxwell's construction 
(dashed line). The construction is justified by analyzing 
Eqs. (fTT|) and at the equal area point where 



exp 



L/3 



dr' i^i - y^{r')) 



= 1. 



(18) 



Here r*^.^ and r1^^^^ denote the value of r at the minima of 
G(/x, r) that correspond to the ordered and homogenous 
phases, respectively. This method, employed for various 
values of /3, yields the first order transition line (thick 
line in Fig.[2b). The right and left dashed lines in Fig.[5]D 
correspond to the stability limit of the homogenous and 
ordered phases, respectively. They define the coexistence 
region where we find three extrema of G {11, r). 

The picture emerging from the derivation above can be 
verified using Monte Carlo simulations, displayed in Fig. 
[21 In the conserving simulation, consisting of processes 
(HI) and ij.{r) is evaluated by measuring the average 
number of ABC triplets in the lattice. The results show 
good agreement with the hydrodynamic solution, thus 
confirming the validity of Eq. (|4]) . In the nonconserving 
simulations, consisting of processes ([l])-([3|), the average 
value of r is measured. The results deviate from the theo- 
retical curve only close to the first order transition shown 
in Fig. [3]d. There, we find different values of r depend- 



ing on whether the simulation was initiated in the fully 
phase separated or homogeneous states. This hysteretic 
behaviour is an indication of a first order transition. The 
agreement with the theoretical results demonstrates the 
applicability of the large L limit considered above to fi- 
nite systems with slow nonconserving processes. 

To conclude, we have studied the generalized ABC 
model with slow nonconserving dynamics. This limit en- 
ables us to derive an exact expression for the LDF of the 
overall density, r, based on the knowledge of the con- 
serving steady-state, p*(a;, r), despite the fact that the 
LDF of pa{x) is not known. In addition, we define the 
chemical potential of the model, which unlike the equi- 
librium one, depends on the details of the nonconserving 
dynamics. Based on this approach we compute the ex- 
act phase diagrams of the conserving and nonconserving 
models. They consist of a second order transition line 
which turns into a first order line at different tricriti- 
cal poiirts in each model. Such ensemble inequivalence 
is typical of equilibrium models with long-range interac- 
tions. This suggests that due to long-range correlations, 
which appear generically in driven diffusive systems 
l3l| . the 'grand-canonical' phase diagram derived follow- 
ing the approach presented above may often differ from 
the corresponding 'canonical' phase diagram. 

The derivation above can be readily applied to other 
driven models that are coupled slowly to an external 
bath. The explicit expression of the LDF of the iroir- 
conserved parameter may, however, depend on the steady 
state properties of the conserving model which are known 
analytically only in a handful of models. We demon- 
strate this approach and the significance of the choice of 
the nonconserving process by considering the generalized 
ABC model with the usual grand-canonical dynamics, 
where the nonconserving process ((3|) is replaced by 



0, 



B 



C 



(19) 



Assuming slow nonconserving dynamics ([S]) and following 
the derivation above yields the same LDF as in Eq. 
but with /i(r) — [log(r) — log(l — r)]//]. As expected, 
different nonconserving dynamics lead to different def- 
initions of the chemical potential. Here, since fi{r) is 
single- valued for any value of /3, the conserving and non- 
conserving model display the same phase diagram. It 
consists of a second order transition line given by Eq. 
dHI) for A 0. 

It would be interesting to study the borderline case of 
7 = 2 and investigate how the picture presented above 
changes when the 'adiabatic' approximation ^ breaks 
down. A similar limit has recently been studied in a 
boundary-driven diffusive model with nonconserving dy- 
namics in the bulk, for which an implicit expression for 
the LDF of the profile was derived [l2j. In a different 
study, an approximate chemical potential has been de- 
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fined for several driven models where the nonconserving 
dynamics is not slow (32Ll33|. 
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